Lattice effects on the spin dynamics in antiferromagnetic molecular rings 
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We investigate spin dynamics in antiferromagnetic (AF) molecular rings at finite temperature in 
the presence of spin-pfionon (s-p) interaction. We derive a general expression for the spin suscep- 
tibility in the weak s-p coupling limit and then we focus on the low-frequency behavior, in order 
to discuss a possible microscopic mechanism for nuclear relaxation in this class of magnetic materi- 
als. To lowest order in a perturbative expansion, we find that the susceptibility takes a Lorentzian 
profile and all spin operators (5*^, ,S^) contribute to spin dynamics at wave vectors g 7^ 0. Spin 
anisotropies and local s-p coupling play a key role in the proposed mechanism. Our results prove that 
small changes in the spatial symmetry of the ring induce qualitative changes in the spin dynamics 
at the nuclear frequency, providing a novel mechanism for nuclear relaxation. Possible experiments 
are proposed. 
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In a crystal of molecular magnets not all physical prop- 
erties and microscopic interactions can be exhaustively 
described by an cfFectivc Hamiltonian containing only 
magnetic terms. The works on magnetization tunnel- 
ing in Mni2 compounds iSi^ clearly demonstrated that 
the localized spins in each molecule interact with lattice 
degrees of freedom, dislocations and organic ligand 
deformations—, changing local field at each magnetic site. 
Spin-lattice interactions revealed thus an essential ingre- 
dient for explaining the relaxation of the macroscopic 
magnetization in ferromagnetic clusters. In the case 
of antiferromagnetic (AF) rings, thermodynamic prop- 
erties such as specific heat and magnetization curves can 
be rather satisfactorily described in terms of a minimal 
iV— spins Hamiltonian containing only magnetic interac- 
tions 
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where periodic boundary conditions are understood. In 
addition to the dominant Heisenberg interaction and to 
the Zeeman coupling, the anisotropy term in Eq. 
(_D ~ — 10~^ J) is an effective description of the local 
field at each magnetic site. This model has been also 
successfully adopted for interpreting the behavior of dy- 
namical observables at the electronic frequency - elec- 
tron paramagnetic resonance (EPR)& or inelastic neutron 
scattering (INS) experiments^. However Hamiltonian ^ 
proves not adequate for explaining nuclear relaxation 
phenomena in AF rings and a coupling with external de- 
grees of freedom is then required: The discrete spectrum 
of the molecular magnet does not allow the relaxation be- 
tween Zeeman-split nuclear levels, due to the mismatch 
between the typical nuclear energy scale (10~^J) and 
the electronic one (ranging from J to gUbH ~ 10^^ J). 
NMR experiments performed on AF rings — suggest the 
relevance of the spin-phonon (s — p) coupling for under- 
standing and explaining the relaxation rate l/Ti 
and. more generally, the spin dynamics at very low fre- 



quency in these systems. In spite of the large amount 
of experimental NMR data available in Feg — and Crg 
— , a comprehensive microscopic description of the relax- 
ation mechanism in AF rings at finite temperature and 
weak magnetic field is still missing. Past investigations 
were mainly focused on the dynamics at high magnetic 
field and zero temperature in connection to the problem 
of tunneling of the Neel vector — . Assuming that the 
Hamiltonian is a good starting effective model for 
the AF rings, a central problem is to clarify which spin 
mechanisms control the low- frequency dynamics, allow- 
ing to observe the Lorentzian profile of the relaxation rate 
1/Ti and its strong frequency dependence^iiS. Moreover, 
it is not fully understood which magnetic property may 
be investigated by the nuclear relaxation measurements, 
due to the mismatch between the nuclear and the elec- 
tronic energy scale. A microscopic description of the s-p 
effects on the spin dynamics should also help the identifi- 
cation of the dccohcrcncc mechanism in this class of ma- 
terials, the interaction of the spins with external degrees 
of freedom strongly limiting possible future technological 
applications—. 

In this Letter we investigate the dynamical spin sus- 
ceptibility in the presence of a spin-phonon coupling in 
order to identify the possible microscopic spin mecha- 
nisms at the basis of the observed dynamics. We derive 
a general expression for the dynamical spin correlation 
functions and then we focus on the very low frequency 
limit, connecting our results to the problem of the nu- 
clear relaxation in AF rings. We will show that, to low- 
est order in the s-p coupling, even small perturbations 
breaking the reflection symmetry of the molecular ring 
have a profound effect on the relaxation rate, due to con- 
tributions coming from the dynamics of all spin operators 
(5^ and 3^,8^) at wave vectors k 7^ 0. In a recent work 
— based on a Master Equation formalism, Santini et al. 
found that in the uj —^ limit, the spin dynamics at low 
temperature, in the presence of a coupling with a ther- 
mal bath, is dominated by the operator at q ~ wave 
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vector. As clearly shown in a related context, however, a 
fully microscopic approach is often essential for clarifying 
quantum effects in mesoscopic spin-phonon systems — . 

Due to s-p coupling, although the nucleus does not 
directly interact with phonons the nuclear relaxation 
rate 1/Ti probes the spin dynamics of the many-body 
system, the relaxation rate being expressed as^ 
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depending on the form of the s-p coupling V{S) and on 
the spin state: a'l = ^^^X=^. Substituting the product 

eigenstates into the formal expression ||2Jl we obtain 
S''^k,io) ^ 
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In Eq. Q A^'^ are geometrical coefficients and S""^{k, uj) 
represents the dynamical spin correlation function at the 
nuclear Larmor frequency ojl ~ 10^^ J. The dynamical 
correlations may be conveniently expressed in Lehmann 
representation as 
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where is a generic single spin operator (fc — i'Kn/N 
n = . . . iV — 1), and \i > and \f > are the many-body 
exact eigenstates of the full Hamiltonian, including spin 
and phonon degrees of freedom, with energy Ei and Ef 
respectively. The system of interest is described by the 
minimal model Hamiltonian 
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where aq are the Bose operators for the (three dimen- 
sional) phonons, ujq the phonon frequency and V{S) is 
a generic spin operator describing the coupling with the 
phonons {s-p coupling constant is implicitly included in 
the definition of the potential V{S)). Assuming that the 
s-p term weakly perturbs the spin system, we express 
the generic state \i > (and |/ >) in Eq. as a prod- 
uct of eigenstates jcr > |m >, |ct > describing the spins 
and |m > the phonon degrees of freedom. To lowest 
order in a perturbative expansion in the s-p coupling, 
we take |cr > as an exact eigenstate of the Hamiltonian 
(^. Then we write the eigenvalue equation for the state 
\i >= \a > \m > and, after projecting onto the spin state 
|(T >, we finally obtain the equation defining the phonon 
state |m >: 

Ea\m > +^^ujqalaq\m > -{- (7) 
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The exact solutions of Eq. (Tj) are products over q of 
shifted harmonic oscillator states, the shift constant 



where \ai > and |cr/ > identify pure spin states with en- 
ergy difference AEif and \n > and |m > are the phonon 
states with phonon occupation number {Uq} and {niq} 
respectively. The key element in Eq. ^ is the overlap 
between the two shifted harmonic oscillator eigenstates 
I < n\m > I, the two shifts depending on the spin states 
|(Ti > and \(7f >. Using an integral representation of the 
delta function in Eq. (jSJ and the explicit expression for 
I < n\m > we obtain 
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where Zp is the phonon partition function and the func- 
tion Iif{oj) is defined (see also—) 
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In the simple situation of a phonon spectrum of the De- 
bye type: ujq = cq for cu < cupi, taking the thermody- 
namic limit and considering the low temperature behav- 
ior Pcup) — > oo, Eq. Hll() becomes 
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where A(u) = J^" dz 



and we have introduced two 
important quantities, the dimensionless frequency ^if = 
\(3{lo — AEif) and the effective coupling 

j,f = 3u;^^<t\V\i>-<f\V\f>]^ (12) 



In the weak s-p coupling limit, 7;/ 
reduces to 
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and Eq. ^ 
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Finally, factorizing the partition function in terms of 
phonon and spin terms Z = Z^Z^ and collecting together 
the previous expressions, we obtain 
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FIG. 1: Dynamical susceptibility S"^{q = vr, oj) as a function 
of the temperature T/J for Crs (Af = 8; 5 = 3/2; J ~ Il.TK) 
calculated according to Eq. 1141 at fixed magnetic field {h = 1 
T,9 = 0.5 rad) and at different frequencies (in units of J). 
Filled symbols refer to the operator, empty to . Inset; 
Exact Diagonalization Results for the spin susceptibility of 
the pure spin system (i.e. no spin-phonon coupling) at lj = 0, 
calculated at two different magnetic fields. 



Eq. (|14|) . valid in the weak coupling regime, reproduces a 
Lorentzian profile and allows to connect in a non-obvious 
way the three different energy scales present in our prob- 
lem: the nuclear, the thermal and the magnetic one. Re- 
gardless the functional form of the spin-phonon interac- 
tion V{S), we first observe that when the initial spin 
state \(Ti > is equal to the final one \af >, the effective 
coupling jif vanishes and lif reduces to a delta func- 
tion at w = 0. Within our formalism, the first non-zero 
contribution at low frequency in Eq. (|14|l is obtained by 
summing over pairs of spin states for which jif =/= 0. If the 
energy gap between the initial and final states is larger 
than the frequency scale (i.e. ^Eif » lu), frequency 
dependence in Eq. H14|) can be neglected, giving thus an 
w— independent relaxation rate 1/Ti. The experimental 
observation of a strong frequency dependence in the 1 /Ti 
suggests that other channels contribute to the relaxation 
mechanism of the nuclear spin. Processes involving states 
originally belonging to the same spin multiplet are not 
good candidates, because the energy splitting created by 
the presence of both anisotropy and field in a generic ori- 
entation is two or three orders of magnitude larger than 
the nuclear energy scale. In order to preserve a frequency 
dependence in Eq. H14|) we arc instead forced to consider 
transitions between pairs of quasi-degenerate spin states 
Eif ^ lol- By examining the spectrum of Hamiltonian 
Q , we found that degenerate states connected by reflec- 
tion symmetry - with momenta q and —q- are the only 
possible candidates for describing the mechanism gov- 
erning the very low frequency dynamics at temperature 
T ^ J, within our approach. However, in order to get 
a non vanishing jif via Eq. H12|) . the s-p coupling op- 



erator V{S) must have a non vanishing matrix element 
between the two states \q > and \ — q >, i.e. the reflec- 
tion symmetry of the ring must be explicitly broken by 
the spin-phonon term. Starting from the Hamiltonian 
the interaction of each single ion with the neigh- 
boring organic atoms provides the required perturbation 
terms. Actually it is well established that ligand groups 
surrounding each magnetic site can assume different, en- 
ergetically equivalent, spatial orientations (quenched dis- 
order), lowering the symmetry of the molecule with re- 
spect to that of a perfect coplanar ringi The positions 
of the atoms determine the strength and the type of in- 
teraction of each spin with the phonons, the spin-phonon 
potential being described in a quite general way by the 
expression 
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The operators 0{Si) describe electrical quadrupole 
interactioniiii, while Ci is the site dependent s-p cou- 
pling which breaks the reflection symmetry of the Hamil- 
tonian A detailed description of the functional form 
of the quadrupolar operators 0{Si) and a quantitative es- 
timate of each local spin-phonon coupling is beyond the 
aims of our work and it is not essential for a qualitative 
understanding of the proposed mechanism for the spin 
dynamics. The main point we want to stress is the exis- 
tence of a site dependent interaction, breaking the trans- 
lation and reflection symmetry of the ring: the degener- 
acy between the two states \q > and | — g > is lifted and 
a non vanishing coupling constant jif (|12|l is generated. 
The nuclear relaxation rate thus acquires a contribution 
from the spin dynamics at wave vector k = 2q showing a 
strong frequency dependence when the splitting created 
by the perturbation is comparable with the nuclear en- 
ergy scale. The differences in the interactions at each 
magnetic site can be usually neglected at energies higher 
than the nuclear one and site dependent parameters do 
not improve the quality of the experimental fiti. On the 
contrary, wc argue that, at the nuclear frequency, very 
small differences in the s-p coupling arc relevant for the 
relaxation rate 1/Ti. We find here a typical aspect of the 
disorder in correlated systems: macroscopic observables, 
as the 1 /Ti , are strongly affected by very small perturba- 
tions, when the perturbation itself becomes comparable 
with the main energy scale of the problem. 

Exact diagonalizations have been performed in order 
to evaluate the low frequency behavior of the dynamic 
correlations in Crg and Fcg molecular rings. In Fig^and 
in Fig. 121 wc report the results for the dynamical sus- 
ceptibility S'^'^{k,Lu) as a function of temperature T/J 
from Eq. for the Crg (fc = tt) and Fee {k = 2n/3) 

compounds. Calculations are performed at fixed mag- 
netic field {h = 1 T, in the direction 9 = 0.5 rad) and for 
different frequencies, including in Ea. H14() pairs of quasi- 
degenerate states. The dynamical susceptibility of both 
and operators shows a maximum as a function of 
the temperature, the value of the maximum decreasing 
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FIG. 2: Same as Fig. [Hfor Fee (iV = 6; S = 5/2; J ~ 28^-). 



and its position moving to higher temperatures when the 
frequency is increased. This behavior quaUtatively repro- 
duces the experimental observations in both materials^. 
The matrix elements in Eq. H14|l are only weakly depen- 
dent on the magnetic field, as can be inferred from the 
results for the spin susceptibility S{k,Lu = 0) of the pure 
spin Hamiltonian (insets of Fig ^ and Fig. |2J). Our re- 
sults support the conclusion that the external magnetic 
field affects the relaxation rate mainly through the Lar- 
mor frequency lul ~ H~. In conclusion we investi- 



gated low-frequency spin dynamics at finite temperature 
in AF rings in the presence of s-p coupling, deriving a 
general expression for the dynamical susceptibility. We 
emphasized the central role of the lattice degrees of free- 
dom, showing that the interaction of the spins with the 
surrounding atoms can generate a contribution to the dy- 
namics at the nuclear frequency from non-diagonal oper- 
ators 5^, q ^ 0. Small changes in the spatial symme- 
try of the ring induce thus novel effects in the spin dy- 
namics at low frequency. A characteristic feature of our 
expressions is that the low frequency dynamics acquires 
a non-zero contribution from all spin operators and 
5*. For isotropic hyperfine coupling A'^'^ oc (5°^'^ (uniform 
magnetic field along direction, radio-frequency field in 
X— direction) l/Ti involves only matrix elements of S^. 
NMR experiments performed on Cr nucleus in Crg or on 
^^Fe in Fcg, should provide an interesting test for detect- 
ing the contribution coming from transverse spin fluctu- 
ations at q ^ 0. Conversely, the proposed mechanism is 
excluded or strongly suppressed when the reflection sym- 
metry is broken at an energy scale larger than the nuclear 
one. By applying hydrostatic pressure it is in principle 
possible to modify the global symmetry of each molecule, 
verifying the role of the reflection symmetry breaking on 
the nuclear relaxation. 
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